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ABSTRACT 

Parameters necessary to analyze the stratified performance of the 
Apollo oxygen tanks include "g" levels, tank elasticity, flow rates 
and pressurized volumes. Methods for estimating "g" levels and flow 
rates from flight plans prior to flight, and from guidance and system 
data for use in the post flight analysis are described. Equilibrium 
thermodynamic equations are developed for the effects of tank 
elasticity and pressurized volumes on the tank pressure response and 
their relative magnitudes are discussed. Correlations of tank 
pressures and heater temperatures from flight data with the results 
of the stratification model developed by C. K. Forester are shown. 
Heater temperatures were also estimated with empirical heat transfer 
equations. Empirical equations were found to yield satisfactory 
agreement with flight data when fluid properties were averaged rather 
than evaluated at the mean film temperature. 


Symbol s 



NOMENCLATURE 


Area 

Tank wall thickness 

Constant in heater temperature sensor lag equation 
Specific heat at constant pressure 
Stefan Boltzman constant 

Constant in Rayleigh number heat transfer equation 

Diameter 

Young's modulus 



Symbols 

F 
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GM 

h 

I 

K 
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M 

MC 

N 

P 

Q 

R 

R a 

r 

T 

t 

U 

V 

a 

3 

e 

P 

* 

0 

a 

u 

a) 


Thrust vector 

Acceleration in Earth gravity units 

Product of gravitational constant and attracting 
body mass 

Enthalpy 

Vehicle inertia matrix 
Thermal conductivity 
Length 
Mass 

Heater thermal mass 
Poly tropic exponent 
Pressure 

Quantity of heat 
Position vector 
Rayleigh number 
Tank radius 
Temperature 
Time 

Internal energy 
Volume 

Vehicle angular acceleration vector 
Coefficient of thermal expansion 
Heater emissivity (0.2 assumed) 

Density 


Thermodynamic property. 


1 

au 

p aP 


ah 

3p 


Thermodynamic property, -p 
Poisson's ratio 
Viscosity 

Vehicle angular velocity vector 
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Subscripts 

b 

eg 

d 

f 

9 

h 

L 

o 

t 

s 


Bulk fluid 
Center of gravity 
Demand 
Thrust 

Attracting body (earth or moon) 

Heater 

Lines 

Reference state or condition 

Tank 

Sensor 


INTRODUCTION AND SUMMARY 

Flight performance of single phase cryogenic tanks (pressure decays 
and heater temperatures) can be reasonably predicted by the techniques 
developed for the Apollo oxygen tanks. Very careful attention is 
necessary in the determination of flight operating conditions 
(accelerations and flow rates) to obtain accurate flight predictions. 

The primary analysis tool used for the tank thermodynamic predictions 
was the numerical model* developed by C. K. Forester for a non- 
isothermal cryogen. It was evident that if this analysis method were 
to provide accurate flight predictions: 

1. Acceleration must be accurately established to model 
the convection processes in the tank, and 

2. Accurate flow rates are required to determine the pressure 
response to heat addition. 

Procedures to calculate these accelerations and flow rates for Apollo 
14 were developed. Early comparisons of the simulation results to 
flight data indicated that additional effects, such as heater thermal 
mass and tank elasticity, needed to be included in the basic solution. 
The basic method of solution did not require modifications. The 
resulting prediction techniques have been correlated to Apollo 14 
flight data and show excellent agreement. 

♦Forester, C. K. , "Pressurized Expulsion of Non- Isothermal Single Phase 
Cryogen", paper presented at the NASA-MSC Cryogenics Symposium, 

May 20-21, 1971. 
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Heater temperatures analyses were also made with an empirical Rayleigh 
number equation for convection. The fluid properties for this simpli- 
fied analysis were averaged instead of evaluated using standard con- 
ventions. This approach resulted in good correlations with flight 
data and provides a convenient tool for heater temperature predictions. 


SYSTEM DESCRIPTION 

The oxygen flow distribution system is shown by Figure 1. The system 
includes check valves which are intended to prevent flow into the 
tanks during normal operation. The isolation valve between tanks 2 
and 3 is normally open and for the Apollo 14 mission was closed only 
during the high flow test. The flow restrictors are capillary tubes 
which produce a pressure drop of approximately 800 psi at 3.5 lbs/hr 
flow rate. The restrictors are the only significant source of 
pressure drop in the system. 

The data available from the system include fluid quantity and pressure, 
and heater temperature for each of the three cryogenic tanks. The 
surge tank is instrumented to provide pressure data only. The flow 
rate to the environmental control system is measured downstream of 
the surge tank and, therefore, includes contributions from all four of 
the tanks. The flow rate to the fuel cells is also measured, but 
can be more accurately determined from the electrical current. 


FLOW RATES 

The total flow from the three tanks during the Apollo 14 mission was 
determined from the fuel cell usage and the flow rate across the 
restrictors to the environmental control system (ECS). The fuel cell 
flow rate was computed using the fuel cell current, because the 
current data are more accurate than the fuel cell flow meters. The 
flow rate across the restrictors during high flow periods was based on 
the restrictor pressure drop calibrations. During low flow periods, 
the restrictor flows were obtained from the ECS flow rate (measured 
downstream of the restrictors) and the net change of mass in the surge 
tank during the period. 

The flow rates from the individual cryogenic tanks were not measured, 
therefore, it was necessary to divide the total system flow among the 
three tanks. The individual tank flow rates were determined from the 
total system flow on the basis of equilibrium tank thermodynamics. The 
pressure differences between tanks were used to determine the check 
valves configuration and to constrain the thermodynamics calculations. 
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FIGURE 1 - APOLLO 14 OXYGEN SYSTEM FLOW SCHEMATIC 





The flow distribution is affected by the heat input to the separate 
tanks, therefore, the heat leaks were estimated from flight data. The 
tank 1 and tank 3 heat leaks were found to be nominal at zero flow 
rate at 90% and 10% quantities, respectively. The tank 2 heat leak 
was not verified, but is believed to have been slightly greater than 
nominally expected. The tank 2 heat leak could not be determined for 
a zero flow condition because the check valve provided to isolate the 
tank leaked. The check valve permitted warm fluid to flow back into 
the tank causing warming of the insulation and increasing the heat 
leak. 


The flow rate distributions were obtained by simultaneous solution of 
the pressure change equations (Appendix) for the tanks supplying the 
system flow. The calculations included the effects of tank elasticity, 
since the elasticity strongly affects the pressure change rate at 
high quantities (Figure 2). The simultaneous solution of equations 
A-19 for two tanks supplying the system demand relates the individual 
tank flow to the total flow. The total flow used for this calculation 
included the flow rate required to pressurize the external line volumes. 
These calculations are simplified if the pressure change rate is known 
and used with the nominal heat leak to determine the individual tank 
flows. The flow rates from this method are in the same ratio as those 
provided by simultaneous solution of the equations. 


The tank 2 check valve leak required special consideration to deter- 
mine the flow into the tank causing pressurization. The flow rate 
into the tank was determined from the volume change required to produce 
the observed pressure change. Using equation A- 1 5 and considering the 
volume change due to a hot bubble as well as tank elasticity, the 
pressure change is: 


dP ♦ 
3t " V 




0) 


where 
density 


dV 

t 


is the volume change due to a bubble of fluid at the line 


dM_ 


During pressurization, the demand flow d is zero. Substituting 


-f L M£) L 


and solving for ^ we have: 
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This method of determining tank 2 inflow during the tank 3 heater cycle 
at 26 hours AET resulted in check valve leakage rates of approximately 
0.05 lbs/hr which is believed to be realistic. 

Thermodynamic analyses were not required to distribute the flows during 
the high flow test, since the isolation valve was closed. The tank 3 
flow rates (Figure 3) during the tests were determined from the restric- 
tor pressure drop while the surge tank valve was open. Flow rates 
were assumed constant while the surge tank valve was closed. The tank 
1 flow rates (Figure 4) were also determined from the restrictor 
pressure drop, but included the fuel cell flow rate when the tank 1 
pressure was greater than the tank 2 pressure. 


TANK ACCELERATION ANALYSIS 


The sources of accelerations in a space vehicle in drifting flight 
include vehicle rotations, thrusts caused by fluid venting, gravity 
gradients and solar pressure. The solar pressure is approximately 
10~ 7 lbs/ft 2 and produces an acceleration of less than 5 x 10"^ "g" 
for the Apollo vehicle. The acceleration due to solar pressure was an 
order of magnitude smaller than the accelerations produced by vehicle 
rotations during a typical Apollo 12 attitude hold period and was, 
therefore, neglected. 

The procedure for the analysis of accelerations during attitude hold 
conditions used rotation rates from guidance data directly for the 
centripetal acceleration. The rotation rates were numerically 
differentiated for the angular acceleration term. The total accelera- 
tion due to rotation is: 

32.174 g = w X W X (R t - R cg ) + a x (R t - R cg ) (3) 

Telemetry data from the digital auto pilot used for the analysis 
includes the three components of the is vector and the calculation is, 
in principle, straightforward. Some difficulty does, however, arise 
due to the angular acceleration. The angular acceleration terms tends 
to dominate the centripetal term because the centripetal acceleration 
depends on the square of the rotation rate. The acceleration term 
also introduces questions of significance due to the short durations 
of application. Typically, the reaction control system jet firings 
cause angular accelerations greater than 2 x 10“4 radians/second 2 , but 
the duration is of the order of 10 milliseconds. This acceleration 
results in a movement of the oxygen tank of about 10“ 7 inches during 
the time the acceleration is applied. This small displacement would 
appear to be negligible, however, the angular accelerations should 
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certainly not be entirely ignored. The approach used was to distribute 
the angular acceleration over time intervals of 10 seconds or greater 
by numerical differentiation of the observed angular rates at the end 
points of the time interval. The time intervals were selected on the 
basis of engineering judgement to adequately characterize the accelera- 
tion events. Although this procedure is arbitrary, the results appear 
to be satisfactory and a better method has not presented itself. Tank 
accelerations during an Apollo 12 attitude hold condition were 
typically 5 x 10"° to 7 x 10“° "g". 

The tank accelerations during stable periods of passive thermal control 
(PTC) were calculated without consideration of angular accelerations. 
During PTC flight modes, the reaction control system is deactivated 
and the vehicle is essentially spin stabilized. For this condition, 
angular accelerations are generally negligible and the centripetal 
acceleration only is significant. 


When the vehicle is in attitude hold in the near vicinity of the 
earth or moon, the gravity gradient acceleration is significant, 
term, which must be added to the rotational accelerations, is: 

2 GM (R„ 


This 


32.174 g = 


R g 4 "V 


V R g 


(4) 


The gravity gradient term is of the order of 10“ 7 "g" for a 100 mile 
altitude earth orbit. Since the magnitude of the gravity gradient 


3 

is proportional to 1 / | R | , the term becomes negligible at distances 

of 2-3 earth radii. The radius vector to the attracting body in the 
vehicle coordinate system is necessary to the gravity gradient calcula- 
tion. This vector can only be determined from the vehicle trajectory 
and inertial platform data. A computer program for the calculation 
of the acceleration including the gravity gradient term derived from 
trajectory data has been developed by NASA-MSC for application to the 
Apollo 14 mission. The average tank acceleration in lunar orbit was 
5 x 10“ 7 "g" for the Apollo 12 mission. 


The tank accelerations for nominal flight conditions can reasonably 
be assumed to be the same as previous flights at similar conditions. 
The high flow oxygen tank tests during the Apollo 14 mission were not, 
however, in nominal flight "g" conditions due to overboard dumping of 
oxygen. The oxygen dumped overboard through a convergent nozzle in 
the command module entry hatch produced a significant thrust and 
vehicle acceleration. 


The thrust of the 6.2 Ibs/hr oxygen dumped overboard was calculated 
as the thrust from a choked convergent nozzle for a stagnation 
temperature of 60°F and a specific heat ratio of 1.4. The resulting 
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thrust of 0.091 pounds is conservative, since no expansion downstream 
of the throat or plume effects were considered. 

The 0.091 pounds thrust produced a linear "g" of 3.65 x 10 -6 for the 
24,985 pounds vehicle weight at the time of the test. The thrust vector 
was not through the vehicle center of mass, therefore, rotational 
accelerations were also produced. The equation for the angular 
accelerations is: 

R f x F = I o ( 5 ) 

and solving for 

a = I" 1 (l? f x T) (6) 

The moments of inertia obtained from pre-flight mass properties data 
were used with equation 6 to predict angular acceleration rates. The 
rotation rates were then calculated from the time required to rotate 
the vehicle through the 5° dead band. Finally, rotational tank acceler- 
ations were calculated from equation 3. The total tank 3 "g" was found 
to be 4.9 x 10 -6 and for tank 1 a "g" of 4.7 x 10" 6 was obtained. These 
accelerations neglected the effects of reaction control system firings 
as well as plume effects and therefore, may be somewhat lower than 
actually experienced. The pre-flight acceleration predictions have 
been used for the high flow test post flight analysis, since accelera- 
tion data reduction has not been completed. 


ANALYSIS CORRELATIONS WITH FLIGHT DATA 

Evaluations of the numerical model and empirical heater temperature 
equations were based on Apollo 14 flight data. The accelerations and 
flow rates required by the analyses were determined using the methods 
described earlier. Comparisons of simulations with Apollo 14 data 
showed good correlation with the flight pressures and temperatures. 

THERMODYNAMIC SIMULATIONS 

Detailed tank thermodynamic simulations were conducted with the 
numerical math model described in the paper* by C. K. Forester. This 
model uses a rectangular (x and y) grid of cells to approximate the 
flow fields in the tank. The simulations were accomplished for a range 
of grid sizes to permit extrapolations of variables to asymptotic limits. 
The effects of tank elasticity and heater temperature sensor laq were 
included. 


*IBID, page 3. 
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The tank pressure and heater temperature results from a simulation of 
one PTC heater cycle are compared with flight data by Figure 5. This 
simulation was made with a heater area of 0.475 ft^ which is the flat 
plate area equivalent to the 0.59 ft^ outer surface of the cylindrical 
heater tube. Since the heater tube is perforated, flow through the tube 
could provide an effective flat plate heater area of 0.95 ft^. Simula- 
tions were conducted for both heater areas to determine which provided 
the most accurate heater temperature. The results of the convergence 
analyses for the areas are shown by Figure 6, The asymptotic limit 
for the heater temperature with the 0.475 ft^ area is within 9°F of the 
flight data. The asymptotic temperature with the larger heater area is 
not in good agreement with flight data. The asymptotic limit of the 
heater on time for the small heater is also in better agreement with 
flight data than the large heater. These results imply that the inside 
of the heater tube is not an effective heat transfer surface. 

The simulation of the tank 3 high flow test is shown by Figure 7. The 
convergence study and the asymptotic limit for the heater temperature 
are shown by Figure 8. The asymptotic heater temperature with the 0.95 
ft^ heater area is 45°F below flight data. The PTC heater cycle 
simulation with this area resulted in a heater temperature 95°F below 
the flight data. Since the 0.475 ft^ heater area provided good agree- 
ment for the PTC heater cycle, the 0.95 ft^ area was expected to 
produce temperatures much lower than flight data. The heater tempera- 
ture simulation results of the PTC heater cycle and the high flow test 
are contradictory. The contradiction may be due to the "g" level which 
was estimated as 4.9 x 10" 7 "g" for the high flow test. A higher "g" 
level during the high flow test could bring the results into agreement. 
The discrepancy could also be due to residual fluid rotations caused 
by the prior PTC period. Resolution of this discrepancy must await 
more accurate determination of the "g" level from flight data. 

EMPIRICAL HEAT TRANSFER EQUATIONS 

Heater temperatures can be determined from the numerical math model, 
but the computer time requirements are excessive for the generation 
of parametric data and routine flight analyses. Empirical heat 
transfer equations were investigated to develop a more convenient 
tool for heater temperature studies. 

The convective heat transfer from a horizontal cylinder is usually 
determined from a Rayleigh number equation. 

|a = n L K AT C ra (R a ) 1/ ‘* (7) 
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FIGURE 5 - PTC HEATER CYCLE SIMULATION 
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HEATER SENSOR TEMPERATURE - °F HEATER ON TIME - MIN 



NUMBER OF X-DIRECTION CELLS 
FIGURE 6 - PTC HEATER CYCLE CONVERGENCE 
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HEATER TEMPERATURE - °F PRESSURE - PSIA 



FIGURE 7 - TANK 3 TEST SIMULATION 


118 




HEATER TEMPERATURE - °F HEATER ON TIME - MIN 
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FIGURE 8 - TANK 3 TEST CONVERGENCE 
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( 8 ) 


The Rayleigh number is determined from: 

D 3 p 2 32.174 9 MT C. 
R a W 


The fluid properties used to evaluate the Rayleigh number are usually 
taken at the mean film temperature. This convention is based on tests 
with simple fluids under 1 "g" conditions. Since the properties of 
supercritical oxygen may vary by an order of magnitude in the boundary 
layer, the properties in the Rayleigh number were averaged instead of 
taken at the mean film temperature. The viscosity, conductivity, and 
density were taken as the average of their values for the bulk 
temperature and the heater temperature. The specific heat was 
evaluated as the difference in the enthalpy at the heater and bulk 
temperatures divided by the temperature difference. The coefficient 
of expansion used was. 


0 = 



(9) 


The radiation from the heater is also significant and was included 
in the complete heat transfer equation. 


$=„LKATC ra (R a ) lA + e C $ (T h 4 - T b 4 )A 


( 10 ) 


Heater temperatures were developed as a function of on time by 
numerical integration of the equation. 



(ID 


where MC is the heater thermal mass of 0.1 BTU/°F. The heater 
temperature sensor lag was included in the integration to provide 
a means of comparison with flight data. The temperature sensor 
response was determined from: 


dT s 

3r 


= c k < T k 


V 


( 12 ) 


The temperature-time results of a simulation of one heater cycle 
using this simplified analysis (Figure 9) are in excellent agreement 
with flight data. Comparisons between actual and predicted heater 
temperatures (Figure 10) for several heater cycles indicates 
predictions within 50°F with the exception of the tank 3 heater cycle 
during the high flow test. As already mentioned, the "g" level 
assumed for this period is questionable. 
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Parametric heater temperatures (Figure 11 ) were generated for easy 
comparisons with flight data. The results of the comparisons made 
so far are promising, however, the full ranges of quantities and 
accelerations have not been investigated. A final determination of 
the empirical constants and the averaging method used for fluid 
properties can not be made, since attitude hold conditions have not 
been investigated. 


CONCLUSIONS 


Conclusions based on correlations with the Apollo 14 high flow tank 
tests and passive thermal control conditions are: 

1. The flight performance of single phase cryogenic tanks 
can be reasonably predicted using the numerical model 
developed. Heat transfer and detailed thermodynamic 
processes (stratification) in the tanks are accurately 
modeled. 

2. Heater temperatures can be predicted from empirical heat 
transfer equations if fluid properties are averaged 
instead of evaluated at the mean film temperature. 

3. The accuracy of performance predictions is dependent 

on the ability to predict or control flight accelerations 
and flow rates. 
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HEATER SENSOR TEMPERATURE 



INDICATED QUANTITY - % 

FIGURE 11 - PARAMETRIC HEATER TEMPERATURES 


APPENDIX 


THE PRESSURE CHANGE EQUATION FOR A 
CRYOGENIC TANK 


The pressure changes in a cryogenic tank resulting from heat addition 
and mass extraction are usually calculated with the assumption that the 
tank is a constant volume container. This assumption causes large 
errors when the fluid is nearly incompressible and the pressure vessel 
is highly stressed. An error in the pressure change calculation is also 
caused by flows not usually measured that are required to pressurize 
plumbing system volumes at ambient temperature. In order to eliminate 
these errors, the pressure change equation for an equilibrium fluid in 
an elastic container has been derived and a method for including the 
external volume effects developed. 

The thermodynamic system is bounded by the inside surface of the 
pressure vessel and is closed at the fluid outlet from the pressure 
vessel. The volume inside the thermodynamic boundary is not constant 
since the pressure vessel is elastic. The outflow velocity is assumed 
to be small enough that the kinetic energy and momentum of the outflow 
are negligible. The conservation equations for mass and energy, 
therefore, determine the system response to heat and mass flows. 

The conservation of mass: 


dM ..do . dV 

3t ■ v at + B 7t 

The conservation of energy: 

at 

Expanding A-2 and using the definition of enthalpy 

Substituting A-l in A-3 and simplifying 

u dU dQ . P dM D dV 

pv ?t vS p^t' p 7t 


(A-l) 


(A-2) 


(A-3) 


(A-4) 
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The internal energy is taken as a function of pressure and density 


dU = 3U dP 9U d£_ 

Bt ~ aF cTt 3p dt 


( A- 5 ) 


Using A-l and noting that M = pV 


dU _ 3U dP . 3U /I dM M dV\ 


Substituting A-6 in A-4 



3U dP 
W (It 


+ 


3U/1 dM M dV 
3p V? clt " 7 (3Ft 



P dM D dV 
p 3t " <Jt 


Solving for ^ and rearranging 



dP .. 

d£ P dM 
3t P 3t 


3U /dM M dV\ 
p 3^Uf " V Si) 


3t 




Rearranging 

dP _ 

3M?( 

P _ 3U> 

P P 3p> 

\ 3U M dV D dV 
1. p 3p” V cTt ~ at 


cTt " 



* „3U 

1>V 3P 

and reducing 

the last 

term 




d( 

dP <T 

1+^(1 
t 3t\p 

3U\ 

' P 3P/+ 

2 dV/3U P \ 

p ^(^- 7 ) 


<it ' 

u3U 

cV 3P 


w3U 

PV 3P 




Now define 


and 



(A-6) 
(A- 7 ) 

(A-8) 

(A-9) 

(A-10) 

(A-l 1 ) 
(A-l 2) 
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(A— 13) 


Substituting A- 1 1 and A- 1 2 in A-10 


dP _ * / dQ . ^ \ 1 dV 

dt • 7 (at + ®3tr f 3t 0 0 * 

Equation A- 13 provides a convenient method for calculating pressure 
response if the volume rate of change is known. The last term is 
zero for a constant volume system. If the container is elastic so 
that the volume change is related to the pressure change, some further 
reduction is possible. 

dP _ /dQ, adM\ 1 dV dP ... 

<Jt V \fit + ^cft / " V cfF cflf pe<t> ( A ‘ 14 ) 


Again solving for ^r 


. 1 

dP V 

It 


1 + 


dM 

t ' "Ht 

r av 

V w p6<t> 


for a spherical tank 


therefore 


1 dV 

VW 


3 r (1-a) 
5bE 


dP 

dt 


1 


(A-15) 


(A- 16) 


(A-17) 


Now the outflow is measured at the end of the distribution lines 
which contain gas at the same pressure as the tank. The flow across 
the thermodynamic boundary must include the flow required to pressurize 
the lines. Two assumptions for determining the flow into the lines 
may be considered: 

1. Fluid expelled from the tanks thermodynamic boundary 
maintains its density while compressing gas in the lines 
either adiabatical ly or isothermally. 

2. Fluid expelled from the tank does not affect the tempera- 
ture distribution from the thermodynamic boundary to 

the system outlet, but the density in the lines changes 
adiabatically or isothermally. 
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The phenomena of assumption 1 can be described by writing the 
poly tropic relationship for the volume of gas in the lines. 



Taking logarithms and differentiating we have: 



J_ dP 
N P 


(A-18) 


(A-19) 


The flow rate into the tank thermodynamic boundary due to the lines 
is therefore: 


dM L 

ir = - 


\ dP 
NP It 


(A-20) 


Where N=1 for isothermal compression and N=1.4 for adiabatic 
compression. 


Now writing A- 14 in terms of the demand and line flow rates. 


dP 

It 




69 dM d 1 dV dP 
V IT ' V IP It 


p 4>e 


Substituting A-20 in A-21 and again solving for 



dP 

It 


we have: 


(A-21) 


(A-22) 


The relationships for assumption 2 are developed by assuming that the 
fluid density in the lines is related to the pressure by the poly- 
tropic exponent. 


p - 



(A- 23) 


Again taking logarithms and differentiating we have: 


d£ _ 1 dP 
p " N P 


(A-24) 
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(A-25) 


Therefore, the flow rate into the tank from the lines is: 

d«i p L Y L dp 

cTt TJP 

dD 

Substituting A-25 into A-14 and solving for -pjr as before we have: 



(A-26) 
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